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ABSTRACT 

Recent radio observations have identified a class of structures, so-called radio 
relics, in clusters of galaxies. The radio emission from these sources is interpreted 
as synchrotron radiation from GeV electrons gyrating in /iG-level magnetic fields. 
Radio relics, located mostly in the outskirts of clusters, seem to associate with 
shock waves, especially those developed during mergers. In fact, they seem to be 
good structures to identify and probe such shocks in intracluster media (ICMs), 
provided we understand the electron acceleration and re-acceleration at those 
shocks. In this paper, we describe time-dependent simulations for diffusive shock 
acceleration at weak shocks that are expected to be found in ICMs. Freshly 
injected as well as pre-existing populations of cosmic-ray (CR) electrons are con- 
sidered, and energy losses via synchrotron and inverse Compton are included. We 
then compare the synchrotron flux and spectral distributions estimated from the 
simulations with those in two well-observed radio relics in CIZA J2242. 8+5301 
and ZwC10008. 8+5215. Considering that the CR electron injection is rather in- 
efficient at weak shocks with Mach number M < a few, the existence of radio 
relics could indicate the pre-existing population of low-energy CR electrons in 
ICMs. The implication of our results on the merger shock scenario of radio relics 
is discussed. 

Subject headings: acceleration of particles — cosmic rays — galaxies: clusters: 
general — shock waves 
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INTRODUCTION 



The presence of energetic nonthermal particles, especially electrons, in clusters of galax- 
ies has been inferred from observations of so-called "radio halos" and "radio relics" (see, 
e.g., ICarilli fc Tavlorl I2OO2I : iGovoni & Ferettil bood ; iFerrari et aDbood : iBriiggen et aDboil 



for reviews). The radio emission from these sources is interpreted as synchrotron radiation 
of cosmic-ray (CR) electrons. The radio halos center roughly in cluster cores and have low 
surface brightness with steep radio spectrum and low polarization. Radio relics, on the con- 
trary, are isolated structures, typically located in the cluster outskirts but within virial radii. 
They often exhibit sharp edges, and most of them show strong polarization. In fact, with 
occasional pairings found in the opposite side of clusters and elongated morphologies, radio 
relics are commonly thought to reveal shock waves in intracluster media (ICMs) produce d 
during mergers (e.g. , lEnfilin et al.lll998l ; iRoettiger. Burns fc Stondll999l ; iMiniati et al.ll200ll ) . 
Unfortunately, relics are found mostly too far from cluster cores for their X-ray signatures to 
be easily detected. So only in a few cases, their association with ICM shocks have been es - 
tablished by X-ray observations (e.g.. lFinoguenov et al.ll2010l ; lAkamatsu &: Kawaharall201l[) . 



More than 40 relics have been identified in radio observations so far ( Nuza et al .1120121 . and 
references therein). Based on the spatial distribution of shocks seen in cluster formation sim- 
ulations, it is predic t ed that coming radio surveys will e asily identify hundreds more (e.g., 



Skillman et al. 



20 111 : IVazza et al.ll20ia iNuza et all 120121 ). 



The observed synchrotron radiation is expected to come from CR electrons with Lorentz 
factors 7 e > 10 4 , spiraling in ~ fiG magnetic fields. The cooling time scale of such CR 
electrons due to synchrotron emission and inverse Compton (IC) scattering does not much 
exceed ~ 10 8 yrs (see equation (J2J)). Advection or diffusion over that time would typically 
be limited to < 100 kpc. So, the electrons have very likely been injected, accelerated or 
re-accelerated close to where they are seen in emission. 

Shocks, believed to be associated to observed radio relics, are obvious candidates for 
the acceleration or re-acceleration of the CR electrons. Suprathermal particles are known 
to be produced as a n inevitable consequence of the formation of collisionless shocks in ten- 



uous plasmas (e.g.. iGarate fc Spitkovskyl |2012| ). If postshock suprathermal particles have 



sufficient rigidity to recross the shock transition, they can be fur ther accel e rated to be- 



come CRs through so- called Diffusive Shock Acceleration (DSA) (jBelll 1 19781 ; iDruryl 11983 
Malkov Druryl l200ll ). Only a very small fraction of inflowing plasma particles are "in- 



jected" from the thermal pool into the CR population. Yet, in strong shocks, a sufficient 
number of CRs reach high energies, so that they extract a substantial fraction of the dissi- 
pated energy, allowing DSA to be efficient. 



Shock waves are indeed common in the intergalactic space (e.g., 



Miniati et al 



2000: 
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Ryu et al.l 120031 ). They are induced by the supersonic flow motions produced during the 
hierarchical formation of the large-scale structure (LSS) in the universe. Those shocks are, 
in fact, the dominant means to dissipate the gravitational energy which is released during 
the LSS formation. They broadly reflect the dynamics of baryonic matter in the LSS of the 
universe, and, indirectly, dark matter. Simulations suggest that while very strong shocks 
form in relatively cooler environments in filaments and outside cluster virial radii, shocks 
produced by mergers and flow motions in hotter ICMs are relatively weak with Mach number 



M < a few (IRvu et al 



Hoeft et al J 1200a 



200 



3; 



Vazza et al 



'frommer et al. 



2006 



2009; iBriiggen et al.l 



2011( 1. 



Kang et al.l 120071 : ISkillman et alJl2008 



At weak shocks, however, DSA should be inefficient. This is expected from the fact that 
the particle energy spectrum associated with DSA is steep when the density compression 
across a shock is small. Also the relative difference between the postshock thermal and flow 
speeds is greater in weaker shocks. Consequently, t he injection from thermal to nonthermal 
particles should be inefficient at weak shocks (e.g.. lKang et al.ll2007l ). At shocks with M < 
a few, much less than ~ 10 -3 of protons are thought to be injected into CRs and much 
less than ~ 1 % of the shock r am pressure be converted into the downstream pressure of 
CR protons ( jKang fc Ryull2010l ). For reference, recent Fermi observations of 7-ray emission 
from galaxy clusters, searching for 7-ray by-products of p — p collisions, limit the pressure 
due to CR protons to le ss than ~ 10 % of the gas thermal pressure there f lAbdo et al. 



2010l ; lDonnert et al.ll2010l ). IACT (Imaging Atmospheric Cherenkov Technique) observations 
c rn x ~ _ " _ " J 1 ' of < 1 — 2 % in core regions of some clusters 



Alaksic et al. 


2012; 


Arlen et al. 


2012) 



Injection and acceleration of electrons are even more problematic at weak shocks. Rel- 
ativistic electrons and protons of the same energy are accelerated the same in DSA, since 
they have the same rigidity. But nonrelativistic electrons of a given energy have substantially 
smaller rigidities than protons, making them much harder to be injected at shocks from the 
thermal pool. As a consequence, the number of electrons injected and accelerated to the 
CR population is likely to be significantly smaller than that of CR protons, and so is the 
pressure of CR electrons at weak shocks. 

Hot ICMs, on the other hand, should have gone first through accretion shocks of high 
Mach numbers aroun d clusters and fi l aments and then through weaker shocks inside those 
nonlinear structures ( IRvu et al.l l2003t iKang et al.l 120071 ). Hence, it is expected that ICMs 
contain some CR populations produced through DSA at the structure formation shocks. 
In ad dition, in ICMs, nonthermal par ticles can be produced via turbulent re-acceleration 
(e.g., iBrunetti fc Lazarianl 120071 . 120111 ). Moreover, secondary CR electrons are also con- 



tinuously generated through p — p collisions of CR protons with thermal protons of ICMs 
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(e.g., iMiniati et al.ll2001r iPfrommer & EnBlinl 12004] ). If radio relics form in media with such 
"pre-existing" CRs, the problem of inefficient injection at weak shocks can be alleviated. 

In this paper, we study DSA of CR electrons at shocks expected to be found in ICMs, 
with and without pre-existing CR electrons. Since the shocks are mostly weak, the CR 
pressure is likely to be a small fraction of the thermal pressure (see iKang fc Ryul 1201 ll ). 
So we apply DSA in the test-particle regime. In the time-asymptotic limit without radia- 
tive losses, the test-particle DSA theory predicts a steady-state distribution of power-law 
for downstream CR electrons, f e ,2(p) °c P q with q = 3cr/(cr — 1), w here a is the density 
compression ratio across a shock, when no pre-existing CR is assumed (jDrurylll983l ). If pre- 
existing CR electrons of a power-law distribution, f e ^ oc p~ s , are assumed, the distribution 
of re-accelerated el ectrons approaches f e p{jp) oc p~ r with r = min(g, s) at large momenta 
(IKang fc Ryull201ll . and also see equation (|UJ)). The power-law distributions of f e ,2(p) trans- 
late into the synchrotron/IC spectra of j„ oc wit h a = (q — 3)/2 or (r — 3)/2 (e.g., 
Zirakashvili fc Aharonian 2007 ; Blasi 201ol ; Kang 2011 ). These properties provide essential 
benchmarks for expected spectral properties 

We perform "time-dependent", DSA simulations of CR electrons for plane-parallel 
shocks, which include the energy losses due to synchrotron and IC processes. Using the 
simulation data, we calculate the synchrotron emission from CR electrons, and model the 
synchrotron flux and spectral distributions from spherical shocks. We then compare the re- 
sulting distributions to t hose of two well-observed radio relics in clusters C IZA J2242. 8+5301 
Ivan Weeren et alJboioh and ZwC10008.8+5215 Ivan Weeren et al.l 120111 ) in details. 



The relic in CIZA J2242. 8+5301 at redshift z = 0.1921 perhaps demonstrates the best 
evidence for DSA at merger shocks. It is located at a distance of ~ 1.5 Mpc from the cluster 
center and spans ~ 2 Mpc in length. The relic shows a spectral index gradient towards the 
cluster center. The spectral index, measured between 2.3 and 0.61 GHz, steepens from —0.6 
to —2.0 across the relic. The relic is strongly polarized at the 50 — 60 % level, indicating 
ordered magnetic fields aligned with the relic. In the opposite, southern part of the cluster, an 
accompanying fainter and smaller relic is found. The relic in ZwC10008. 8+5215 at z = 0.1032 
is found at a distance of ~ 0.85 Mpc from the cluster center and has a linear extension of 
~ 1.4 Mpc. It also shows the steepening of the spectral index towards the cluster center. 
The spectral index, measured between 1382 and 241 MHz, changes from —1.2 to —2.0 across 
the relic. The polarization fraction is less with < 25 %. It also has an accompanying relic of 
a linear extension of ~ 290 kpc in the opposite, western side of the cluster. 

In Section 2 we describe the numerical method and the models for magnetic field, dif- 
fusion, electron injection, and pre-existing CR electron population. We present analytic 
evaluations for some features in the CR electron energy spectrum and synchrotron emission 
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spectrum in Section 3. The results of simulations are presented and compared with obser- 
vations of the previously mentioned radio relics in Section 4. Summary follows in Section 
5. 



2. DSA SIMULATIONS OF CR ELECTRONS 



2.1. Numerical Method 



We simulate DSA of CR electrons at gasdynamical shocks in one-dimensional plane- 
parallel geometry. Shocks in ICMs, especially merger shock s, are expected to pe rsist over 



> 10 9 yrs, a substantial fraction of the cluster lifetime (e.g., ISkillman et al.ll201l[ ). On the 
other hand, the time scales over which electrons are accelerated and cool are much shorter 
(see Eq. [2] below). So we assume that the shock structure remains steady. Assuming 
that the CR feedback to the flow is negligible at weak shocks in the test-particle limit, the 
background flow, u, is given by the usual shock jump condition. Then, the time-dependent 
evolution of the CR electron distribution, f e (t, x,p), which is averaged over pitch angles, can 
be followed by the following diffusion convection equation, 



dg e + u d£k _ Idu f dg e 
dt dx 3 dx \ dy 



d_ 

dx 



<x,y) — 



d (b e 
+ PTT- —9e 

ay \p z 



where g e = p A f e , y = hi(p/m e c), m P is the electr on mass, c is the speed of light, and n(x,y) 
is the spatial diffusion coefficient ( jSkillingj Il975h . 

Here, b e (p) = (4e A /9mlc 6 )B^ s p 2 represents the cooling of CR electrons due synchrotron 
and IC losses in cgs units, where e is the electron charge. The "effective" magnetic field 
strength, £?| ff = B 2 + -E>cbr> includes the equivalent strength of the cosmic background 
radiation with -Bcbr = 3.24 /xG(l + z) 2 at redshift z. The cooling time scale for electrons is 
given as 



trad(7e 



P 



be(p) 



9.8 x 10 7 yrs 



5 fiG 



7e 

10 4 



(2) 



where 7 e is the Lorentz factor of CR electrons. 



The equation in (ED) is solved using t he test-particle version of the CRASH (Cosmic-Ray 
Amr SHock) code (see lKang et al.ll201lL for details). 
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2.2. Models for Magnetic Field and Diffusion 



Here, shocks are assumed to be gasdynamical for simplicity, that is, magnetic fields do 
not modify th e background flow of the shock. In ICMs, magn etic fields of an inferred strength 
of order /zG (jCarilli fc Taylor! l2002r iGovoni fc Ferettil 12004 ) are dynamically un important, 



since their energy density is less than ~ 10 % of the thermal energy density (e.g.. iRyu et al. 
20081 ). However, magnetic fields, especially in the downstream region, are the key that 
governs DSA and the synchrotron cooling and emission of CR electrons. Theoretical studies 
have shown that efficient magnetic field amplification via resonant and non- resonant wave- 
parti cle interactions is an integral part of DSA at strong shocks (ILucek fc Belli 120001 iBell 
2004 ). In addition, magn e tic fields can be amplified by turbulent motions behind shocks 
( IGiacalone fc Jokipiil 120071 ; llnoue et al.l 120091 ) . Yet, these plasma processes are complex and 
their roles are not yet entirely certain, especially at weak shocks. So here we adopt a simple 
model in which the magnetic field strength is amplified by a constant factor of \ across the 
shock, that is, B 2 = xB\- Hereafter, we use the subscripts '1', and '2' to label conditions in 
the preshock and postshock regions, respectively. 

For k, we adopt a Bohm-like diffusion coefficient with weaker non-relativistic momentum 
dependence, 



K{x,p) 



P 



m P c 



(3) 



where k* = m e c 3 /(3eBi) = 1.7 x 10 19 cm 2 s l (B\/l fiG) 1 in the preshock region and 
K 2 = K \/x m the postshock region. 



2.3. Injection of Electrons 



As pointed in Introduction, the injection of electrons is expected to be much harder than 
that of protons in the so-called thermal leakage injection model. Because complex plasma 
interactions among CRs, waves, and the underlying gas flow are not fully understood, it is not 
yet possible to predict from first princ i ples how particles are injecte d into the first-order Fermi 
process (e.g.. iMalkov fc Druryl l200ll ; iGarate fc Spitkovskyl |2012| ). In addition, postshock 
thermal electrons, which have gyro-radii smaller than those of thermal protons, need to be 
pre-accelerated to several times the peak momentum of thermal protons, p Pi th, before they can 
re-cross the shock transition layer. Here, p Pjt h = y/2m p fcBT 2 , where T 2 is the postshock gas 
temperature and ks is the Boltzmann constant. Recently several authors have suggested pre- 
acceleration mechanisms based on plasma interactions wit h fluctuating magnetic fields that 
are locally quasi-perpendicular to the shock surface (e.g. iBurgessI l2006t lAmano fc Hoshino 
20091 ; iGuo fc Giacalondl2010l ; iRiquelme fc Spitkovskyl 1201 ll ). But the detailed picture of the 
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electron injection is not well constrained by plasma physics. Observationally, the ratio of CR 
electron number to proton number, K e / p ~ 0.01, is commonly inferred for strong supernova 
remnant shocks, since about 1% of the Galactic CR flux near a GeV is due to electrons 



( [Reynolds! 120081 ) . But this ratio is rather uncertain for weak shocks under consideration. 



So here we adopt a simple model in which the postshock electrons above a certain 
injection momentum, p in j = QinjPp.th, are assumed to be injected to the CR population. 
Here, is a parameter that depends on the shock Mach number and t urbulent magnetic 
field amplitude in the thermal leakage injection model (IKang fc Ryull2010l ). The CR electron 
number density or, equivalently, the distribution function at p in j at the shock location x s , 
/ e (x s ,pi n j), will be scaled to match the observed flux of radio relics (see Sections 3.1 and 
4.2). 



2.4. Pre-existing CR Electrons 

We consider the population of pre-existing CR electrons, along with that of freshly 
injected electrons at the shock. However, the nature of pre-existing CR electrons in ICMs 
is not well constrained. If they were generated at previous, external and internal shocks, a 
spectral slope of s ~ 4.0 — 5.3 is expected for M > 2, close to the acceleration site. However, 
since their lifetime in equation ([2]) is much shorter than that of host clusters, it is unlikely that 
they are directly responsible for the pre-existing electron population we consider here. Any 
pre-existing CR electron should be locally produced, possibly either through p — p collisions 
of CR protons with thermal protons or via turbulent re-acceleration of some populations 
(possi bly including p — p secondary electrons), as noted in Introduction. iPetrosian fc East 



( 120081 ) have shown that turbulent injection of CR electrons from the thermal pool in ICMs is 
unlikely. The sl ope of protons re-accelerated by turbulence should be cl ose to s ~ 4 (see, e.g., 
Chandran] 2005 ). but that of electrons is strongly modified by coolings ( Brunetti &: Lazarian 



2007 



20111 ) . The slope of secondary elect rons from p — p collisions would be roughly s 



4/3(s p — 1) (IMannheim fc Schlickeiserlll994l ). where s p is the slope of CR protons, so typically, 
s ~ 4 — 6. In summary, pre-existing CR electrons may contain many different populations 
with different degrees of radiative cooling and may not be represented by a single power-law. 

For simplicity, here we adopt a power-law form, 



f e ,i(p) = f x 



pre 



V 



(4) 



with slope s, as the model spectrum for pre-existing CR electrons. In modeling of specific 
radio relics, the value of s will be chosen as s = 2a Q bs + 3, where a Q b s is the observed mean 
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spectral index. The amplitude, / pre , is set by the ratio of upstream CR electron pressure to 
gas pressure, Ri = PcRe,i/P g ,i- Here, Ri is a parameter that will be scaled to match the 
observed fluxes of radio relics (see Sections 3.1 and 4.2). 



3. ANALYTIC EVALUATIONS 

We first consider some features in the CR electron energy spectrum and synchrotron 
emission spectrum for plane-parallel shocks, to provide analytic estimations for the simula- 
tion results presented in the next section. 



3.1. Basic Features in CR Electron Spectrum 



In the test-particle regime of DSA, the distribution of freshly injected and accelerated 
electrons at the "shock location" can be approximated, once it reaches equilibrium, by a 
power-law spectrum with super-exponential cutoff, 



feM ~ fi 



p v 



cxp 



Peq 



(5) 



where q = 3a/ (a — 1) (jKangll201ll ). In the case that B 2 = <jB>i, that is, the jump in the 
magnetic field strength across the shock is assumed to be same as the density jump, \ = a, 
and K2 = Ki/cr, the cutoff momentum, which represents the balance between DSA and the 
radiative cooling, becomes 

1/2 



Peq 



1 1 

rate u* 



B 1 



The corresponding Lorentz factor for typical merger shock parameters is then 

1/2 

, / ?/,„ \ / Hi \ 

7 e , eq w 2 x 10 y q 



(6) 



-1/2 



3000kms-V \B*+B* 



(7) 



Hereafter, the magnetic field strength is given in units of /zG. The acceleration time for 
electrons to reach p eq , so the time for the equilibrium to be achieved, is estimated as 



t eq « (2.4 x 10 4 yrs) q^B^B^ + B^)^ (, 



,3000 km s" 1 . 

This is much shorter than the typical time scale of merger shocks, > 10 9 yrs. For t > t eq , the 
DSA gains balance the radiat ive losses a nd the electron spectrum near the shock location 
asymptotes to a steady-state (jKang)l201ll ). 
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With pre-existing CR electrons given in equation (j4|), the electrons distribution at the 
shock location can be written as the sum of the pre-existing/re-accelerated and freshly in- 
jected/accelerated populations, 



feM 



(«-«) 



1 - 



P 

Pinj 



—q+s 



fv Te \ Pinj J + 



P 

Pinj 



exp 



when s ^ q 
when s = q. 



(9) 



(IKang fc Ryull201ll ). The relative importance of pre-existing to freshly injected populations 
depends on / pre and / in j, as well as on the slopes s and q in our model. For the sake 
of convenience, hereafter we will use the term "injected" electrons for those injected at 
the shock and then accelerated by DSA and the term "re-accelerated" electrons for those 
accelerated from the pre-existing population. 

We here define the CR electron number fraction, £ e = ricR e ,2/ n e,2, as the ratio of CR 
electron number to thermal electron number in the postshock region. Here ncR e ,2 includes 
CR electrons accelerated from both the pre-existing and freshly injected populations. Con- 
sidering that the C R proton number fraction is likely to be £ p < 10~ 4 at weak shocks 



( IKang &: Ryul 120101 ) and K e / p ~ 0.01, £ e ~ 10~ 6 could be regarded as a canonical value. 
We note that the resulting radio emission is linearly scaled with both £ e and the preshock 
gas density, ni, in the test-particle regime, so the combined parameter, ni£ e , can be treated 
as a free parameter. We here fix the preshock gas density, n\ = 10 -4 cm -3 , as a fiducial 
parameter, but vary £ e to match the observed fluxes of radio relics. Another measure is the 
ratio of postshock CR electron pressure to gas pressure, R2 = -PcRe,2/-P g ,2, which depends on 
both £ e and the slopes q and s. In modeling of specific radio relics in Section 4.2, we will 
determine the set of values for £ e , R2 and R\, that matches the observed level of radio flux. 

If we ignore for the moment the modest influence of continued DSA downstream of 
the shock, we can follow the electron population that advects downstream by solving the 
following equation : 

rln f)n 

0, (10) 



dt dy 



where d/dt = d/dt + ud/dx and V = —b e (p)/p = —Ce y . Here, C = (4e /9m 4 c 6 )5g ff is a 
constant. This is basically the equation for downward advetion in the space of y = ln(p/ra e c) 
due to radiative cooling. The general solution of the equation is 

P 



g e (p,t)=G(e-v -Ct) = G 



11 



. 1 _ t/ t ra( j , 

where t ra d = 1/Ce y is the electron cooling time scale. This provides the approximate distri- 
bution of CR electrons at the distance d = 11$ downstream from the shock, where 112 is the 
downstream flow speed. 
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For instance, if the distribution function of the "injected" electrons at the shock location 
(d = 0) is the power- law spectrum, g e (p, 0) = gmj(p/Pmi)~ q+i , the downstream spectrum can 
be approximated as 

(12) 

It should be straightforward to apply the same approximation to the full spectrum given 
in equation flU]). In Figures 1 and 2, we compare the distributions described by equation 
(|TT|) with those from time-dependent DSA simulations, demonstrating that equation (ITT]) 
provides reasonable approximations to the solutions of full DSA simulations (see Table 1 for 
specific model parameters). 



9e{p,d) = g hi ] 



P 



d/u 2 t vaxi )p 



3.2. Basic Features in Synchrotron Emission Spectrum 

Since the synchrotron emission from mono-energetic electrons with 7 e has a broad peak 
around z/ pea k ~ 0.3(3e5 /47im e c)^, for a given observation frequency, z/ Q b s , the greatest 
contribution comes from electrons of the Lorentz factor, 

7 ,^ Kl . 26x ^(-^-) 1/2 (_^)" 1/2 ( 1 + z) v 2 (13) 

Using equations (J2J) and (fl3l) . the cooling time of the electrons emitting at z/ Q b s can be 
estimated approximately as 

t-«8.7 X 10-,»(^)(^ ; )- V, (l + ,)-V». (14) 

The cooling length behind the shock, w 2 t ra d, then becomes 

^-Mio^) (C) (^r /2(i+ * ri/2 <i5) 

Note that B^ S 2 /B^ 2 ~ 15 — 25 for the model parameters considered here. Again, t rad is 
shorter than the typical time scale of merger shocks, > 10 9 yrs. So L rad should represent the 
width of radio emitting region at ^ b s behind plane-parallel shocks. In radio relics, however, 
the observed width is constrained by both L ra( j and the projection angle of spherical shocks 
(see Section 4.2). 

The cutoff energy in the electron spectrum due to the radiative cooling decreases linearly 
with the distance from the shock location, that is, 7e iCU t oc e? -1 , as expected from equation 
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([2]) and shown in Figure 1. At the farthest downstream point, d = u 2 t, where t is the shock 
age, the cutoff energy becomes 



7e , br (t) « 9.82 x 10 2 



t 



10 9 yrs 



5 fiG 



(16) 



If the electron distribution function at the shock location has a power-law form, n e (x s , j e ) cx 
then the volume-integrated electron spectrum downstream steepens by the power-law 



index of one, i.e., N £j2 oc 7, 



-(r+1) 



for 7 e > 7e,br- It is because the width of the s patial 



distri bution of electrons with 7 e decreases as 7 e 1 ( jZirakashvili fc Aharonianl 12007 ; 



Kang 



20111 ). As a consequence, the "volume-integrated" synchrotron spectrum from aged electrons 



has a spectral break, i.e., an increase of the spectral index a by +0.5, at 

-2 



3 eB 2 2 

z/ br = 0.3- 7 br 

47r rripC 



6.1 x 10 6 Hz 



t 



10 9 yrs 



5 fxG J \5 /iG 



B, 



cff,2 



(17) 



So the shock age may be estimated from the break frequency u^, if the magnetic field 
strength is known. 



4. RESULTS OF DSA SIMULATIONS 

4.1. Plane-Parallel Shocks 

The model parameters of our simulations for plane-parallel shocks are summarized in 
Table 1. Here, z is the redshift, c S) i is the preshock sound speed, M is the shock Mach number, 
u 2 is the postshock flow speed, s is the power-law slope of pre-existing CR electrons, and B 2 
is the postshock magnetic field strength. The model name in the first column includes the 
values of M, B 2 , and s; for models without pre-existing CRs, "I" (injection only) is specified. 
For instance, M4.5B7I stands for the model with M = 4.5, B 2 = 7 /iG, and injected CR 
electrons only (no pre-existing CRs), while M2B2.3S4.2 stands for the model with M = 2.0, 
B 2 = 2.3 /iG, and s = 4.2. For the preshock magnetic field strength, B\ = 1 /x G is adopted for 



all m odels, which is close to the typical quoted value in cluster outskirts (see iBriiggen et al. 



201 ll . and references therein). Once Bi < Bqbr, the IC cooling dominates, and the exact 
value of B\ is not important in our models. The model parameters are chosen to match the 
observed properties of radio relics in clusters CIZA J2242. 8+5301 and ZwCl 0008.8+5215 
(see the next subsection for details). For example, M = 4.5 or s = 4.2 is chosen to match 
the observed spectral index, a = 0.6, of the relic in CIZA J2242. 8+5301, and M = 2 or 
s = 5.4 is chosen to match a = 1.2 of the relic in ZwCl 0008.8+5215. For reference, the 
shock compression ratio is a = 3.5 for M = 4.5 and a = 2.3 for M = 2. The values of u 2 and 
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B<2 are chosen to match the observed width of the relics, since they determine the cooling 
length as shown in equation ffl5|) . 

Figure 1 shows the CR electron distribution at different locations downstream of the 
shock, after it has reached the steady state, for M4.5B3.5I, M2B7S4.2, M2B2.3I and M2B2.3S5.4 
models. For the comparison of different models, here the postshock CR electron number frac- 
tion is set to be £ e = 10~ 6 , which sets the vertical amplitude. The injection-only models 
exhibit the power-law distributions with cutoffs due to the cooling, as discussed in the pre- 
vious section. In M2B7S4.2 model, the electrons accelerated from the injected population 
are important only at low energies (7 e < 10 2 5 ) and they dominate in terms of particle num- 
ber, because the "injected" spectrum is much softer than the "re-accelerated" spectrum (i.e. 
q> s). The electrons accelerated from the pre-existing population, on the other hand, dom- 
inate at higher energies including 7 e ~ 10 4 and they are most relevant for the synchrotron 
emission at v ~ 1 GHz. The slope of the accelerated spectrum at high energies is similar to 
that of the pre-existing spectrum, which is consistent with equation (Q. On the contrary, in 
M2B2.3S5.4 model with s ~ q, the "injected" electrons are negligible even at low energies. 
This difference comes about, because with similar numbers of pre-existing CRs, the ampli- 
tude / pre is larger in M2B2.3S5.4 (with s = 5.4) than in M2B7S4.2 (with s = 4.2). The 
numbers of injected electrons should be similar in the two models, because the shock Mach 
number is the same. Note that the re-accelerated spectrum flattens by a factor of ln(p), as 
shown in equation ([9]), because s ~ q in this model. 

The left column of Figure 2 shows the spatial profile of the electron distribution func- 
tion, g e (7 e ,x), at two specific energies (j e ) as a function of the downstream distance for 
the M4.5B7I, M4.5B3.5I and M2B7S4.2 models. For each model the Lorentz factors are 
calculated for z/ obs = 0.61 GHz and 2.3 GHZ according to equation ({TBI . The upper/lower 
curves represent g e of the lower/higher values of 7 e , respectively. The right column of Figure 
2 shows the synchrotron emission, j u (x), at u ohs = 0.61 GHz (upper curves) and 2.3 GHZ 
(lower curves). The solid lines show g e and j u calculated from the DSA simulation results, 
while the dashed lines show the approximate solutions calculated with equation ( 1TT1) . The 
figure demonstrates that the lower energy elections advect further from the shock before 
cooling than higher energy electrons, so the lower-frequency radio emission has larger widths 
than the higher-frequency one. According to equation ffTo^) . the cooling lengths of the elec- 
trons emitting at 0.61 and 2.3 GHz are L ra d ~ 40 and 20 kpc, respectively, in the three 
models. 
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4.2. Modeling of Radio Relics 

As noted above, it should be sufficient to employ the plane shock approximation to 
compute the distributions of CR electrons and their emissivities as a function of the distance 
from the shock surface. In observed radio relics, however, radio emitting shells are likely 
to be curved with finite curvatures along the observer's line of sight (LoS) as well as in 
the plane of the sky. So in modeling of radio relics, the curved shell needs to be projected 
onto the plane of the sky. In that case LoS's from the observer will transect a range of 
shock displacements, and this needs to be taken into account whe n computing the observed 



brigh tness distribution of model relics. Following the approach of Ivan Weeren et al.l ( 12010 



201 ll ). we consider a piece of a spherical shell with outer radius R s , subtended along the LoS 
from to — ip so for the total angle of 2ip. Then, R s and the projection angle ip are the 
parameters that fix the shape of the curved shell to be projected onto the plane of the sky. 
The synchrotron emissivity, j u (erg cm -3 s _1 Hz -1 str -1 ), at each point behind the curved 
shock is approximated as that downstream of plane-parallel shocks discussed in the previous 
subsection. Since we do not consider the polarization of synchrotron emissions here, so, for 
simplicity, the magnetic field lines are assumed to lie in the plane of the sky; that is, the 
angle between the magnetic field vectors and the LoS is fixed at 90°. 

The synchrotron intensity is calculated by integrating the emissivity along the LoS, 
Iv{ r ) — J judl (erg cm -2 s _1 Hz -1 str" 1 ), where r is the distance behind the projected shock 
edge in the plane of the sky. The bound of the path length, I, for given r is determined by 
R s and ip. Then, the observed flux is estimated, assuming a Gaussian beam with e- width, 
9, as 

S u (r) ^I u (r)n6 2 (l + z)- 3 , (18) 

where v = z/ obs (l + z). 

Figure 3 shows the profiles of the synchrotron flux, S u (r), at z/ Q b s = 0.61 GHz (left 
column) and the spectral index, a = —dlnS u /dlnv, estimated with the fluxes at v ohs = 0.61 
and 1.4 GHz (right column) for the M4.5B7I, M4.5B3.5I and M2B7S4.2 models, which are 
designed to reproduce the radio relic in CIZA J2242. 8+5301. The flux is calculated with the 
beam of 6 2 = 0i#2/(4 In 2), 6182 = 16.7" x 12.7 ". They are compared w ith the "deconvolved" 



profile of observed flux taken form Figure 4 of Ivan Weeren et al.l (120101 ) (filled circles). Since 



the observed flux is given in an arbitrary unit in their paper, we scale it so that the peak 
value of S v (r) becomes 5 mJy, which is close to the observed value (private communication 
with R. J. van Weeren). The radius of the spherical shock is set to be R s = 1.5 Mpc and 
two values of projection angle, ip = 10° and 20°, are considered. The observed profile is 
well fitted by the three models, if i/j — 10° is taken. In M4.5B7I and M4.5B3.5I, different 
values of 112 are assumed to match the observed width (see Table 1). The observed value 
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of the spectral index at r = 0, a = 0.6, is reproduced either in the injection-only models 
with M = 4.5 or in the model with pre-existing CRs with the slope s = 4.2, as noted in the 
previous subsection. 

For the fiducial preshock particle density of n% = 10~ 4 cm -3 , the values of the postshock 
electron CR number fraction required to match the peak flux of 5 mJy are £ e = 7.6 x 
1(T 8 , 2.3 x 1(T 7 , and 2.6 x 10~ 7 for M4.5B7I, M4.5B3.5I, and M2B7S4.2, respectively. In 
M2B7S4.2 model the ratio of the pressure of pre-existing CR electrons to gas pressure far 
upstream is Ri ~ 6.7 x 1CT 5 . Those values of £ e and R\ are modest enough that they probably 
are not in conflict with the values expected in clusters. Our results demonstrate that if the 
pre-existing electron population is considered, the radio relic in CIZA J2242. 8+5301 can be 
reproduced even with weak shocks of M ~ 2 or so. We note that R\ is a model parameter 
that sets the amplitude, / pre , of the upstream population, while the fraction £ e is the outcome 
of DSA of both pre-existing and injected electrons. As noted in Figure 1, in M2B7S4.2 model 
the fraction £ e is determined mostly by the "injected" population at low energies, while the 
radio emission is regulated mostly by the "re-accelerated" population at 7 e ~ 10 4 . So we 
should obtain the similar radio flux even with a much lower injection rate for this model, 
and the resulting £ e could be much smaller than the current value of 2.6 x 10~ 7 . 

We point that the radio relic in CIZA J2242. 8+5301 is subtended in the plane of the 
sky over the angle of ~ 60 — 70°. This means that the surface of the shock responsible for 
the relic should be highly elongated with the aspect ratio of ~ (60 — 7O°)/(2 , 0) = ~ 3 — 3.5 
when ip = 10° is adopted. It would not be trivial, if not impossible, for such structure to 
be induced in merger events in clusters. Or the relic may actually consist of a number of 
substructures, which is hinted by the variations in the observed flux profile along the arc in 
the plane of the sky. 

The left column of Figure 4 shows the synchrotron flux profiles at u ^ s = 1.38 GHz, 
for M2B2.3I and M2B2.3S5.4 models, which are designed to reproduce the radio relic in 
ZwC1008.8+5215. The flux is calculated with 9 2 = 0i0 2 /(41n2), Q X Q 2 = 23.5" x 17.0". We 
note that this beam size is fine enough that the convolved profiles with a Gaussian beam 
(dotted and long-dashed lines) are very similar to the unconvolved profiles (solid lines ) . The 
profiles are compared with the observed profile given in Figure 16 of lvan Weeren et al.l (1201 if ) 
(filled circles). Again the observed flux is given in an arbitrary unit, so it is scaled at 5 mJy 
at the peak. The right column shows the profiles of a, estimated with fl uxes at */ n h s = 0-24 



and 1 .38 GHz, along with the observed a also taken from Figure 16 of Ivan Weeren et al. 



( 120111 ) (filled circles). The shock radius is assumed to be R s = 1.0 Mpc and two values of 
projection angle, ip = 25° and 30°, are considered. The two models shown are the same 
except the existence of pre-existing CR electrons in M2B2.3S5.4 model. In M2B2.3S5.4, the 
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"re-accelerated" population dominates over the "injected" population. Yet, the two models 
give similar profiles of S u and a. We see that in our m o dels ip = 30° gives good fits to the 
observed profiles of S u and a, while Ivan Weeren et al.l (120111 ) argued that ip = 22° seems 
to give a reasonable fit. Note that they adopted u 2 = 750 km s _1 and B 2 = 2 fiG, giving 
-^rad = 40 kpc, while in our models u 2 = 1100 km s _1 and B 2 = 2.4 fiG, giving L ra d = 57 
kpc. 

For the assumed value of ni = 10 -4 cm -3 , the postshock CR electron number fraction 
required to match the peak flux of 5 mJy is £ e = 2.1 x 10~ 4 for M2B2.3I, which is six times 
larger than £ e = 3.3 x 10~ 5 for M2B2.3S5.4. This is because the spectral shapes of CR 
electron spectrum below 7 e < 10 2,5 are different in the two models (see Fig. 1 and discussion 
in the previous subsection). The number fraction of CR electrons for M2B2.3I seems too 
large, consi dering that the pos tshock proton CR number fraction is likely to be £ p < 10~ 4 
for M = 2 (IKang fc Ryull2010l ). In M2B2.3S5.4, on the other hand, the ratio of upstream 
CR electrons pressure to gas pressure is R\ ~ 1.2 x 10~ 3 . This seems to be marginal, that 
is, not inconsistent with expected values, considering that the ratio of CR proton pressure 
to gas pressure is < 10~ 2 — 10 _1 in ICMs as noted in Introduction. But we should point 
that the values of £ e and Ri in these two models are dominated by low-energy CR electrons 
with 7 e < 10 3 (see Figure 1), which do not contribute much to the synchrotron radiation 
observed in radio relics. So if the "injected" population in M2B2.3I consists of electrons 
with 7 e > 10 3 only, the required values of £ e could be reduced by a factor of ~ 10, easing 
down the constraint. Since we do not understand fully the plasma interactions involved in 
the pre-acceleration and injection of electrons at the shock, the detailed spectral shape of 
those low energy electrons are very uncertain. 

The top panels of Figure 5 show the profiles of the intensity, I v {r) = J j v dl, at 6 cm 
(5 GHz), 20 cm (1.5 GHz), and 91 cm (0.33 GHz) in arbitrary units as a function of r 
for the M4.5B3.5I, M2B7S4.2 and M2B2.3I models. Here, the projection angle is set to 
be ip = 30°. Since the emissivity j u decreases downstream of the shock, while the path 
length increases with r, the profiles of I u exhibit non-monotonic behaviors. For example, 
the profiles at 6 cm show a slightly concave turnover before it decreases abruptly at r w 200 
kpc. The middle panels show the spectral indices, a| (solid lines) calculated between 6 
and 20 cm and (dashed lines) calculated between 20 and 91 cm, when the projection 
angle is set to be ijj = 10°, 20° and 30°. The general trend is the increase of a 20 and 
as we move away from the projected shock edge at r = 0, reflecting the effects of radiative 
cooling. Also a 20 > af(, that is, the slope is steeper at higher frequencies. The bottom 
panels show the color-color diagram of a 2 ^ versus a| . The rightmost point ( 
corresponds to the projected shock edge. Away from the edge, the loci move towards the 
lower left direction. In both middle and bottom panels, the spectral slopes also show a slightly 
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concav e turnover for large projection angles of ip = 20° and 30°. Recently, Ivan Weeren et al. 



( 120121 ) reported the color-color diagram for the so-called "Toothbrush" relic in cluster 1RXS 
J0603. 3+4214, which shows a spectral behavior that is consistent with the cooled electron 
population downstream of the shock. 



5. SUMMARY 

In an effort to refine our understandings of radio relics in clusters of galaxies, we have 
performed time-dependent, DSA simulations of CR electrons and calculated the synchrotron 
emission from CR electrons for plane-parallel shocks. The energy losses due to synchrotron 
and IC have been explicitly included. Weak shocks expected to be found in ICMs have been 
considered. Both the cases with and without pre-existing CR electrons have been considered. 
The relevant physics of DSA and cooling is well represented by plane-parallel shocks, since 
the time scales over which electrons are accelerated and cool are much shorter than the 
lifetime of merger shocks in clusters and the radio emission is confined to a region of small 
width behind the shock front. We then have modeled the synchrotron flux and spectral 
distributions from spherical shocks by approximating them with plane-parallel shocks and 
projecting to the plane of the sky for the angle from +ip to — ip along the LoS. For the 
specific models which are designed to reproduce radio relics in clusters CIZA J2242. 8+5301 
and ZwC10008. 8+5215, we have compared the resulting distributions with observed ones in 
details. 

The main results are summarized as follows: 

1) The CR electron spectrum becomes steady, after the DSA gains balance the radiative 
losses. The spectrum at the shock location is well approximated by a distribution with super- 
exponential cutoff at p eq , fe,2(p) oc exp (— p 2 / pl q ) ■ The full expressions of f e ,2(p) and p eq are 
given in equations (jUJ) and (JB}. 

2) The spectrum of the downstream CR electrons that have cooled for the advection 
time, t = d/u2, can be approximated with g e (p,d) = G\p/(1 — d/u2t r£L d)} at the distance 
d from the shock location. Here, G is the functional form of the spectrum at the shock 
location of d = 0. The synchrotron emission from this analytic formula provides a reasonable 
approximation to that calculated using DSA simulation results (see Figure 2). 

3) Both the models of M — 4.5 shock without pre-existing CR electrons and M = 2 
shock with pre-existing CR electrons of _f e ,i p~ A ' 2 may explain the observed properties 
of the radio relic in CIZA J2242. 8+5301. The postshock electron CR number fraction of 
£ e ~ 10~ 7 in the injection-only model or the ratio of upstream CR electrons pressure to gas 
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pressure of R± ~ several x 10~ 5 in the model with pre-existing CRs are required to explain 
the observed radio flux of several mJy. Those values of £ e and Ri are modest enough to 
be accommodated in typical clusters. But the surface of the shock responsible for the relic 
should be highly elongated with the aspect ratio of ~ 3 — 3.5. It would not be trivial for 
such structure to be induced in merger events in clusters. 

4) The radio relic in ZwC10008. 8+5215 may be explained by the models of M = 2 shock 
with or without pre-existing CR electrons. However, in the injection-only model, £ e > 10 -4 , 
required to explain the observed radio flux of several mJy, is probably too large for the weak 
shock of M — 2. On the other hand, in the model with pre-existing CRs, R\ ~ 10~ 3 , required 
to explain the observed flux, seems to be marginal, that is, not inconsistent with expected 
values in clusters. In the model, then, the origin of such pre-existing electron population is 
an important topic, but beyond the scope of the present paper. 

5) The color-color diagram of vs has been presented behind the projected shock 
edge. It includes an important information about the evolutionary properties of the post- 
shock electrons. Due to the effect of the projection with limited subtended angle along the 
LoS for spherical shocks, the diagram behaves differently for different projection angles. So it 
may provide an independent way to estimate the projection angle, which is a key parameter 
in modeling of radio relics. 
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Table 1. Parameters for Plane-Parallel Shock Simulations 



Model Name 


z 


c 8 ,i 


M 


u 2 


s 


B 2 


Cluster 






(km s _1 ) 




(km s _1 ) 




{p.G) 




M4.5B7I 


0.1921 


7.8 x 10 2 


4.5 


1.0 x 10 3 




7.0 


CIZA J2242.8+5301 


M4.5B3.5I 


0.1921 


6.0 x 10 2 


4.5 


7.7 x 10 2 




3.5 


CIZA J2242.8+5301 


M2B2.3S4.2 


0.1921 


1.25 x 10 3 


2.0 


1.1 x 10 3 


4.2 


2.3 


CIZA J2242.8+5301 


M2B7S4.2 


0.1921 


1.25 x 10 3 


2.0 


1.1 x 10 3 


4.2 


7.0 


CIZA J2242.8+5301 


M2B2.3I 


0.103 


1.25 x 10 3 


2.0 


1.1 x 10 3 




2.3 


ZwC10008.8+5215 


M2B2.3S5.4 


0.103 


1.25 x 10 3 


2.0 


1.1 x 10 3 


5.4 


2.3 


ZwC10008.8+5215 
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log(7 e ) l0g(7 e ) 

Fig. 1. — The distribution function, g e = f e> 2{le)lti °f CR electrons at 5 different locations 
downstream of the shock for the M4.5B3.5I, M2B7S4.2, M2B2.3I, and M2B2.3S5.4 models 
(see Table 1 for model parameters). In this figure the postshock CR electron number fraction 
is set to be £ e = 1CT 6 for comparison of the models. In all the models, solid lines are the 
results of DSA simulations. In M4.5B3.5I and M2B2.3I, dashed lines show the approximate 
solution given in equation ffTTl) . In M2B7S4.2 and M2B2.3S5.4, the distribution of the 
pre-existing population of CR electrons is shown with dotted lines. The five downstream 
locations are: d = 0.30, 1.4, 5.8, 29, 86 kpc in M4.5B3.5I, d = 0.34, 1.7, 6.1, 34, 100 kpc in 
M2B7S4.2, d = 0.42, 2.0, 8.1, 40, 120 kpc in M2B2.3I and M2B2.3S4.2. 
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Fig. 2. — Left: The distribution function, g e (x), of CR electrons with 7 e 's given in equation 
( TT3"j) for z/ Q bs = 0.61 and 2.3 GHz in the downstream region for the M4.5B7I, M4.5B3.5I, and 
M2B7S4.2 models. Right: The synchrotron emissivity, j u (x), at z/ Q b s = 0.61 and 2.3 GHz 
for the same three models. Solid lines are the results of DSA simulations, and dashed lines 
show the approximate solution given in equation ( TTTj) and the emissivity calculated using 
the solution. 
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Fig. 3. — The synchrotron flux, S u , at 0.61GHz and the spectral index, a, between 0.61 GHz 
and 1.4 GHz for the M4.5B7I, M4.5B3.5I, and M2B7S4.2 models. Spherical shocks with 
radius R s = 1.5 Mpc are assumed, and two projection angles, ip — 10° (s olid lines) and 20° 



(dashe d lines), are considered. Filled circles are the data points taken from Ivan Weeren et al. 



(120101 ) for the radio relic in CIZA J2242. 8+5301. The flux is scaled so that the peak has 5 
mJy. The required values of the postshock electron CR number fraction, £ e , and the ratio of 
upstream and downstream CR electrons pressure to gas pressure, R\ and R 2 , are shown. 
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Fig. 4. — The synchrotron flux, S u , at 1.38GHz and the spectral index, a, between 0.24 
GHz and 1.38 GHz for the M2B2.3I and M2B2.3S5.4 models. Spherical shocks with radius 
R s = 1.0 Mpc are assumed, and two projection angles, ip — 25° (solid lines) and 30° (dashed 
lines), are considered. Dotted and long-dashed lines are the results convolv ed with a Gausian 



beam with the e- width of 12 . Filled circles are the data points taken from Ivan Weeren et al. 



(120111 ) for the radio relic in ZwC10008. 8+5215. The flux is scaled so that the peak has 5 
mJy. The required values of the postshock electron CR number fraction, £ e , and the ratio of 
upstream and downstream CR electrons pressure to gas pressure, R\ and R2, are shown. 
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Fig. 5. — Top: The intensities at 6 cm (dashed lines), 20 cm (dotted), and 91 cm (solid) 
in the case of if; = 30° for the M4.5B3.5I, M2B7S4.2 and M2B24I models. Middle: The 
spectral indices a\ (solid lines) calculated between 6 cm and 20 cm, and afl (dashed lines) 
between 20 cm and 91 cm for three different values of the projection angle, ip. Bottom: The 
color-color diagrams of ct^ vs a^o f° r three different values of i/j. 



